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A modification of the cavity technique for axion dark matter detection is proposed in which the 
cavity is driven with input power instead of being permeated by a static magnetic field. A small but 
detectable fraction of the input power is pumped by the axion field to a receiving mode of frequency 
uji when the resonance condition ui\ = ojo ± rn a is satisfied, where ojo is the frequency of the input 
mode and m a the axion mass. The proposed technique is found to provide an attractive approach 
to dark matter axion detection, especially for high axion masses. 

PACS numbers: 95.35.+d 



Shortly after the Standard Model of elementaryjiar- 
ticles was established, the axion was postulated [1] to 
explain why the strong interactions conserve the discrete 
symmetries P and CP. Further motivation for the exis- 
tence of such a particle came from the realization that 
cold axions are abundantly produced during the QCD 
phase transition in the early universe and that they may 
constitute the dark matter Moreover, it has been 
claimed recently that axions are the dark matter Q be- 
cause cold axions predict a specific phase space distri- 
bution for the halos of isolated disk galaxies which is 
supported by observation. These considerations provide 
strong motivation to search for axion dark matter. 

Axion properties depend mainly on a single parameter 
f a , called the axion decay constant. In particular the 
axion mass (h = c = 1) 
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and its coupling to two photons: 
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Here a{x) is the axion field, E(x) and B(x) the electric 
and magnetic fields, a the fine structure constant, and 
5 7 a model-dependent coefficient of order one. Although 
f a was first thought to be of order the electroweak scale, 
its value is arbitrary a-priori [4J. However, the limits 
from unsuccessful axion searches in particle and nuclear 
physics experiments combined with constraints from stel- 
lar evolution imply f a > 3 • 10 9 GeV Furthermore, 
the requirement that cold axions do not overdose the uni- 
verse implies f a < 10 12 GeV [2.]. The constraints leave 
open an 'axion window': 2 • 10 -3 < m a ^ 6 ■ 1CP 6 eV. 

It was pointed out in ref. [6] that dark matter axions 
may be detected on Earth through their conversion to mi- 
crowave photons in a cavity permeated by a strong mag- 
netic field. Using this principle and building on earlier 
efforts 0, the AD MX experiment has reached sufficient 
sensitivity to detect halo axions at the expected density 
[3,0]. When SQUIDs are installed as front-end ampli- 
fiers and its cavity is cooled to 50 mK, ADMX will 



be able to detect, or set limits on, halo axions at even a 
fraction of the halo density. The present detector would 
find halo axions immediately (in a minute or so) if the 
axion mass were known. Unfortunately the axion mass 
is poorly constrained. The value of the axion mass for 
which axions provide the observed density of cold dark 
matter (Q a = ^cdm = 0.23) is only known to be of order 
10 -5 eV with large uncertainties, at least a factor 10 in 
either direction. 

The cavity technique is most convenient in the 500 
MHz to 1 GHz range (corresponding to 2.068 < m a < 
4.136 fJ,eV) because those are typical resonant frequencies 
for the lowest TM mode of cavities large enough (fraction 
of one (meter) 3 ) to produce a signal, given that available 
(large volume) magnetic field strengths are of order 8 T. 
Only the lowest TM mode of a cylindrical cavity with 
axis in the direction of the magnetic field has a sizeable 
coupling to halo axions. To extend the search to higher 
frequencies one must fill the available magnetic field vol- 
ume with an increasing number of smaller cavities, kept 
in tune and in phase during the search. Although the 
principle of this method has been demonstrated [llj . it 
presents a challenging engineering problem in practice. 

A possible variation on the cavity permeated by a 
static magnetic field technique is a cavity driven with 
power Po at the frequency loq of one of the cavity modes, 
hereafter called the input mode. Let uj\ be the frequency 
of another cavity mode, called the signal mode. When 
the resonance condition lo\ — loq ± E a is satisfied, where 
E a is the energy of halo axions, the steady state power 
Pi in the signal mode through the coupling of Eq. ([2]) is 
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where g = Qa and Qi are the quality factors of 

the cavity in the input and signal modes, and p a is the 
fraction of the local axion dark matter density with en- 
ergy dispersion 5E a < 5<jJ\ = <jJ\/Q\. The input power is 
assumed to be monochromatic. The dimcnsionlcss form 
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factor C is 
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where e is the dielectic constant inside the cavity. C 
is invariant under the interchange of input and signal 
modes. Eq. ([3]) implies therefore that, all other things 
being equal, it is best to drive the cavity in the lowest 
frequency mode of any given pair. 

Let v and 5v be respectively the velocity and veloc- 
ity dispersion of the axions. Eq. (j3|) assumes that the 
de Broglie wavelength (m a w) _1 is large compared to the 
cavity's linear dimensions. Furthermore, the energy dis- 
persion must be less than the bandwidth of the signal 
mode. The axion energy is the sum of rest mass and ki- 



netic energy: E a = m a (l 
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). Hence SE a ~ m a v8v. 



It is often assumed that the velocity distribution of halo 
axions is Maxwellian with velocity dispersion 5v ~ 300 



km/s, in which case SE a 
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However an isolated 



disk galaxy such as the Milky Way continually accretes 
the dark matter surrounding it and the resulting flows 
do not thermalize over the age of the universe 12]. The 
infall flows produce peaks in the velocity spectrum at any 
given spatial location. Furthermore, the flows form caus- 
tics. Caustics are surfaces in physical space where the 
dark matter density is very large. Evidence was found 
for caustic rings in the Milky Way and in other isolated 
disk galaxies. The evidence is summarized in ref. [l3j |. 
Because an observer in a galactic halo may be close to a 
caustic, his local dark matter velocity distribution may 
be dominated by a single flow or a single pair of flows. 
The evidence for caustic rings in the Milky Way implies 
that our local halo density is in fact dominated by a single 
flow, dubbed the 'Big Flow', with known velocity vector, 
density of order 1 GeV/cc, and velocity dispersion less 
than 50 m/s [Hj]. The energy dispersion of the Big Flow 
is thus less than approximately 10~ 10 m a . 

Superconducting radio frequency cavities have quality 
factors as large as 10 10 [15j. Consider however that, even 
if there are axion flows with negligible velocity disper- 
sion, increasing the quality factor indefinitely does not 
generally yield the best search strategy because high Q 
implies a long time to excite the signal mode. The time 
necessary for the power in the signal mode to rise to 90% 
of the steady state power Pi is tQ ~ &± where V\ = ^ . 
The total time spent at a given cavity tune is t = tQ + td 
where td is the measurement integration time. The sig- 
nal to noise ratio of the search is given by the radiometer 
equation [l6j 
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where T n is the total noise temperature (thermal noise 
from the cavity plus electronic noise from the receiver 



chain) and B is the bandwidth. In searching for cold 
flows, the best strategy is to make the bandwidth as nar- 
row as possible, i.e. B = l/t d - Eq. (O implies then that 
td oc 1/Pi cx 1/Qi- The optimal Qi, the one that max- 
imizes s/n for given t, is such that tQ = td and hence 
t = Therefore, to allow a search over a factor two 

v\ 

in axion mass per year, Q\ should not be larger than 

i 

approximately 10 8 ((fjjz) 2 ■ 

When B = l/td, signal averaging is not possible and 
the noise is exponentially distributed, as in the HiRes 
channel of the ADMX experiment Q. To have less than 
100 false positives due to noise fluctuations per factor of 
two in frequency range covered, the probability that such 
a fluctuation exceed the signal must be less than approx- 
imately 10 2 /Qi per cavity bandwidth. For Qi ~ 10 8 and 
exponentially distributed noise, this requires s/n ~ 15. 
Eq. ([5]) with B = t^ 1 = v\/Q\ implies the minimum 
detectable signal power 



Pi = 2.1 • 10~ 20 W 
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To determine how much input power P will be needed, 
it is useful to write 
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<7 7 = 0.36 in the DFSZ axion model, whereas g 7 = —0.97 
in the KSVZ model. Combining Eqs. ©, ® and Q, 
one obtains 
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for the benchmark values of T n , s/n, g y and p a used 
above. Note however that a large fraction of the power 
Po can be recycled and need not contribute to the energy 
cost. Indeed the power losses in the cavity are the sum 
of losses in the walls and losses through the hole through 
which the signal is extracted: Qq 1 = Q^n + Qhoic- ^ 
for example Qwaii = 10 10 and Qhoio = 10 s , the power 
dissipated in the cavity is only one hundred of the in- 
put power circulating through the cavity. The remaining 
99% may in principle be separated from the signal and 
recycled. Assuming that this is done, we have: 
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with all benchmark factors explicitly included. The en- 
ergy stored in the cavity is of order 0.45 MJ. The field 
strengths must be less than 20 to 25 MV/m for the qual- 
ity factor to be as high as 10 10 . The cavity volume must 
then be larger than approximately (6 meter) 3 . 
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Consider a cylindrical cavity, i.e. with constant cross- 
section relative to an axis z. The cavity modes are labeled 



as usual as TE mnp , TM mnp and TEM p [ljj where p refers 



to the z-dependence of the mode. Using Eq. one finds 
C = unless either the input or the signal mode, or both, 
is TE. For TEM p <-» TE m¥ , the form factor is 
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if p + p' is odd. ip defines the TEM mode through E = 
— V<^(x, y) sin \ and <j> defines the TE mode through 

Vcj)(x, y) sin ■ L is the length of the 

= _pL and w' = 



cavity, 5 is its cross-sectional area, w = -^fij- 



For TM mnp o TE m / n / p /, the form 




FIG. 1: Cross-section of a superconducting cavity designed to 
search for axion stimulated TEM — ► TE transitions. 
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if p + p' is odd. -0 defines the TM mode through 



B = ip^xV#r,2/) cos(^). w = 7 2 + (^)^ 
and a;' = ^ 
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Finally, for TE m „„ <H- 
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if p + p' is odd. In all three cases, C = is p + p' is 
even. Generally speaking, C is largest when the numbers 
characterizing the input and signal modes differ as little 
as possible. In particular, setting p' = p + 1 maximizes 
C with respect to p' for given p. 

The cavity may be tuned by introducing a variable 
amount of helium into it. Since the dielectric constant 
of liquid 4 He is 1.05, it is possible to tune the cavity by 
2.47% in this way. The remaining tuning is achieved by 
changing p. Indeed, on resonance 
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Changing p by one unit while keeping pi 
changes m a by 
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Provided L > 41-^-, a gapless axion mass range is covered 
by changing p (coarse tuning) and by introducing a vari- 
able amount of helium in the cavity (fine tuning). There 



is no difficulty with mode crossings since the helium shifts 
all the resonant frequencies proportionately. For given 70 
and 71, the allowed tuning range is m a!max ~ 71 — 70 to 
TCia, min — 77 However, m a .min is achieved only when 
ujo » m a which implies a low Pi/Po ratio; see Eq. ([3]). 
To keep a high P\/Pq ratio, one searches only over a 
range of a factor 2 or 3 with one cavity. To cover the 
next higher axion mass range, 71 is increased by physi- 
cally modifying the cavity. 

Let us present a specific cavity design for TEM — > TE 
conversion. TM «-> TE and TE «-> TE offer additional 
options to be explored in the future. Fig. 1 shows the 
cross-sectional shape of a cavity which is partitioned into 
communicating layers, with short dimension d, in such a 
way that there is a single series (p = 1, 2, 3, .. ) of TEM 
modes. Let us call TEoi p ' the set of TE modes for which 
7 = ^ m the d << b limit; see Fig. 1. The form factor 
for TEMp «-» TE i,p+i is in that limit 
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where lu (a/) is the angular frequency of the TEM (TE) 
mode. Combining Eqs. © and (TT5")) . one finds the total 
input energy required to search the range < m a < 7: 



E ~ 150 kW year /(X) 
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where I(X) = jX 4 — X 2 + § — In X. For example, since 



7(2) = 1.44, the energy cost to cover a factor 2 in axion 
mass is 215 kW year for the benchmark factors indicated 
in Eq. (ITBl . Note that there is no frequency dependence 
on the RHS of that equation other than that implicit in 
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T n and Q wa ii- To cover the range from 0.5 to 1 GHz, 
one builds a cavity with d = 15 cm. For the cavity to 
be tunable by helium and choice of p, it is sufficient that 
L > 41d = 6.15 m. To cover the range from 1 to 2 GHz, 
d is reduced by half, and so forth. 

Although motivated by the expectation that part of the 
local dark matter is in cold flows, the proposed technique 
is useful even if there is no structure at all in the axion 
energy distribution. One may search for a signal of width 
SE a ~ 10 _6 m a by adding of order 10~ 6 Qi successive 
bins of the output spectrum. Since the number of co- 
added bins is of order 100, the signal to noise is reduced 
by a factor of order 10. This loss in sensitivity is partly 
compensated by the fact that the co-added noise is nearly 
Gaussian distributed as a consequence of the central limit 
theorem. A search with s/n ~ 5 is suitable then, whereas 
s/n r~j 15 is needed in the narrowband search. Also the 
broadband search is sensitive to all of the local axion 
density whereas the narrowband search is sensitive to 
cold flows only. Any remaining loss in sensitivity may be 
compensated for by increasing the input power. Searches 
with different assumptions on the energy distribution of 
dark matter axions can be carried out simultaneously, as 
in the ADMX experiment. 

To extract a signal of a few times 10~ 20 W from a 
cavity driven with a few times 10 7 W of input power, 
albeit at a different frequency, may be a challenging task. 
The input power is not exactly monochromatic since its 
frequency is swept in time. However, filters to suppress 
the input power at the frequency of the signal mode may 
be installed. At the output of the cavity the transmitted 
input power must be separated from the signal since it 
would otherwise overwhelm the first stage amplifier in the 
receiver chain. This can be done by installing between 
the detector cavity output and the first stage amplifier 
a circulator followed by a filter cavity (two in series may 
be needed) with a resonant frequency at w\ but none at 
ujq. The filter cavity transmits the signal while reflecting 
the input power toward the circulator whence it becomes 
available for recycling. 

False positives will be produced by fluctuations in the 
noise and by fake signals. We already discussed noise 
fluctuations. One source of fake signals is electromag- 
netic power that leaks into the cavity from the environ- 
ment. This problem is the same as in ADMX where 
such false positives are called 'environmental peaks' and 
where the problem is known to be manageable. If need 
be, environmental peaks may be eliminated altogether 
by placing the whole detector in a Faraday cage. False 
positives may also be produced by anything, other than 
the axion field, that pumps power from the input mode 
to the signal mode. Non-linearities in the cavity response 
may do this when ujx/ujq is a ratio of integers. However, 
because the quality factor is very high, the integers in 
the ratio are normally very large and the effects of non- 
linearities are suppressed then. Furthermore, the ratio 



lji/uiq is constant during tuning by helium. It ought to 
be possible to avoid special w\/u)q ratios which cause a 
problem. Cavity vibrations may pump input power to 
the signal mode. However this effect ought to be small 
because there is normally negligible power in vibrations 
at GHz frequencies. 

We conclude that, assuming a number of technical 
challenges can be overcome, the proposed technique pro- 
vides an attractive approach to axion dark matter detec- 
tion, especially for high axion masses. 
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